Finite Temperature Transitions in Large Magnetic Field in Dipolar Spin Ice 
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We use Monte Carlo simulations to identify the mechanism that allows for phase transitions in 
dipolar spin ice to occur and survive for applied magnetic field, H, much larger in strength than 
that of the spin— spin interactions. In the most generic and highest symmetry case, the spins on one 
out of four sublattices of the pyrochlore decouple from the total local exchange+dipolar+applied 
field. In the special case where H is aligned perfectly along the [110] crystallographic direction, spin 
chains perpendicular to H show a transition to q = X long range order, which proceeds via a one to 
three dimensional crossover. We propose that these transitions are relevant to the origin of specific 
heat features observed in powder samples of the Dy2Ti2 07 spin ice material for H above 1 Tesla. 
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In 1997, Harris et al. found the surprising experimen- 
tal result that, despite overall ferromagnetic interactions 
of approximately 2 K in strengh, the Ho2Ti207 magnetic 
pyrochlore material fails to order down to a temperature 
of 50 mK pj. In Ho2Ti207 the Ho 3+ magnetic moments 
reside on a lattice of corner-shared tetrahedra (Fig. ^) . It 
was argued that the strong local Ising anisotropy forcing 
the Ho 3+ moments to point along the local (111) cubic 
directions frustrates the ferromagnetic interactions 0,0- 
The number of ground states in this frustrated Ising sys- 
tem is macroscopic, giving a residual zero temperature 
entropy close to that estimated by Pauling for the proton 
disordered state in common hexagonal water ice lh [3; 
hence the name spin zee |. Specific heat measurements 
on Dy2Ti20j Ja, another Ising pyrochlore material, and 
Ho2Ti 2 07 [fabi have provided strong evidence that the 
low temperature state of these systems indeed possess an 
entropy close to that found by Pauling for water ice. 

For water I/j, applied hydrostatic pressure introduces 
constraints that lift the proton degeneracy, giving a rich 
pressure-temperature phase diagram and a multitude of 
crystalline ice phases ||. This leads naturally to the 
question of how to introduce constraints in magnetic spin 
ice and to search for analogous ordering behavior. Harris 
et al. []J, and Ramirez et al. [f| first explored this ques- 
tion by investigating the behavior of Ho2Ti207 Q and 
Dy2Ti207 0] under applied magnetic fields. In partic- 
ular, specific heat measurements on powder samples of 
Dy2Ti207 revealed three prominent features [5j. For a 
field H larger than 1 Tesla (T), a sharp peak suggesting 
a phase transition occurs at a field-independent temper- 
ature of 0.35 K, surviving up to the largest field consid- 
ered (6 T). Another rounded peak at 1.2 K first appears 
at 0.75 T and also remains observable up to 6 T, al- 
though it is less sharp than the 0.35 K feature. Finally, 
a third peak at 0.5 K first appears at 1 T, where it is 
quite sharp, but disappears for H > 3 T. Over the past 
six years, much theoretical effort has been devoted to the 
search for a microscopic explanation of the field-induced 
spin correlations responsible for these specific heat fea- 




FIG. 1: Portion of the pyrochlore lattice, with the [110] and 
[112] cubic crystallographic axes. The double headed arrow 
shows a [111] spin on a site of the a sublattice that is decou- 
pled from the shown H (solid arrow). Sites on the a FCC 
sublattice (solid black circles) are shown. 

tures in powder samples of Dy2Ti207 0. 

Building on the work of Harris et al. on the effect of 
a [110] field (see Fig. [TJ on Ho 2 Ti 2 7 Q, Ramirez et 
al. suggested that some of these features might be due 
to phase transitions arising in crystallites that accidently 
happen to be oriented such that some of the Ising spins 
in the unit cell are perpendicular to H. These "field- 
decoupled" spins would then be free to interact among 
themselves, giving transitions at temperatures that are 
field-independent up to a very large field 5]. A similar 
proposal had been made earlier for the case of garnet sys- 
tems 0. Experimental evidence 01 suggests that the 0.5 
K feature in powder samples arises due to a multicritical 
point at a temperature of 0.5 K and field of 1 T along the 
[111] crystallographic direction, which is related to the in- 
teresting phenomenology of kagome ice |ll| . However, to 
date, there is no concensus as to whether or not the spe- 
cific heat features at 0.35 K and 1.2 K are caused by the 
magnetic field directed on crystallites of particular orien- 
tation 0,E3,0|. In this paper we address this question 
by using the dipolar spin ice model, where Ising spins 
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on the pyrochlore lattice interact via nearest-neighbor 
exchange and long range magnetic dipole— dipole inter- 
actions ^3- We find a generic scenario for field- induced 
phase transitions which survive even when the field is 
much larger than the combined exchange and dipolar in- 
teractions. 

It is now well established that the nearest-neighbor 
exchange in Ho 2 Ti 2 7 6] and Dy 2 Ti 2 7 |3 El E3 
is antiferromagnetic and that the spin-ice phenomenon 
originates from long-range magnetic dipole-dipole inter- 
actions rather than from nearest-neighbor ferromagnetic 
exchange as originally proposed 0,0- The magnetic mo- 
ment, /i e ff, of Ho 3+ and Dy 3+ in the above materials is 
/ieflt ~ 10/iB, giving a dipolar coupling constant, D, at 
nearest- neighbor distance of approximately 1.4 K j(ill ,'>l|. 
The pyrochlore lattice is a non-Bravais lattice with four 
atoms per unit cell (e.g. a white tetrahedron in Fig. 1) 
with each sublattice making a regular FCC lattice. The 
dipolar spin ice model Hamiltonian is |l5| : 
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The first term is the Zeeman interaction between the 
spins and the magnetic field H. The second and third 
terms are the exchange and long-range dipolar inter- 
actions, respectively. The spin vector Sf — afz a at 
R" = Hi + r a labels the Ising moment of magnitude 
|S"| = 1 at FCC unit cell coordinate site Ri and tetrahe- 
dral sub-lattice site coordinate r a (a = 1,2,3,4), where 
the local < 111 > Ising axis is denoted by z a and the Ising 
variable is erf = ±1. The vector R? b = R? - Rf con- 
nects spins Sf and S b . If one restricts exchange coupling 
to nearest-neighbor spins, Ji, a ;j,b = J\ represents the ex- 
change energy and D the dipolar energy scale. Ji, a ;j,b > 
is ferromagnetic, Ji. a -j,b < is antiferromagnetic (Ji > 
and D = in the nearest-neighbor spin ice model 0,0)- 
Previous comparison between Monte Carlo simulations 
of H with experimental specific heat data has found that 
Ji « -3.72 K for Dy 2 Ti 2 7 [3 and J x w -1.65 K for 
Ho 2 Ti 2 7 0, with D w 1.41 K for both materials. In 
addition, we assume strictly classical spins and neglect 
quantum fluctuations induced by the applied transverse 
field which admix the excited crystal field states with the 
ground state doublet. This is a reasonable approximation 
for fields smaller than 10 T, since the energy gap to the 
lowest excited states is ~ 300 K in Dy 2 Ti 2 7 |18|. 

Consider first the problem of spins in the dipolar spin 
ice model, subject to an external magnetic field H, from 
the perspective of symmetry. For H / 0, the simplest 
scenario for a phase transition arises when one of the 
four FCC sublattices, call it a (see Fig. 1), is decoupled 
from H. In other words, the field is perpendicular to 
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FIG. 2: Specific heat from MC simulations with only the (a) 
FCC sublattice occupied by spins. The values J = —3.72 K 
and D = 1.41 K appropriate for Dy2Ti20 7 yield a transition 
at 0.65 K. Adding a small J3 = —0.022 K shifts this transition 
to a temperature that agrees with experiments (see text). L 
is the linear size of the system, the total number of spins in a 
pyrochlore lattice is iV = 16L 3 . Inset: Simulations of the full 
pyrochlore lattice with H(f9) (see text). The system behaves 
like the a FCC lattice in the limit of large (Inf.) field. 



the Ising quantization direction z a such that H • z a = 0. 
With the condition H-z a = 0, there exists one completely 
free rotational degree of freedom in the orientation of the 
crystal with respect to H, which allows all spins belong- 
ing to the a sublattice to interact among themselves and 
participate in a phase transition. One obvious field orien- 
tation that satisfies this condition is the [112] direction, 
first investigated in Ref. [l9|. Measurements on a single 
crystal of Dy 2 Ti 2 7 with H aligned along the [112] direc- 
tion revealed two broad specific heat peaks and no sign of 
a phase transition ^H- For H larger than 0.5 T or so, an 
anomaly develops at T-J°^ ~ 1.7 K, independent of field 
strength, while a higher temperature anomaly occurs at 
•^112 ex I H I . We confirmed the results of Higashinaka et 
al. via Monte Carlo simulations on the spin ice model 
for Dy 2 Ti 2 7 . In particular, using the above J\ and D 
values for Dy 2 Ti 2 7 , we find that for an infinite [112] 
field a broad Shottky-like specific heat anomaly occurs 
at ~ 1.7 K, with no sign of a sharp phase transition. 

We now explain the results in Ref. ^3- ^ n the limit 
of strong [112] field, the spins on the three a ^ a sub- 
lattices that have a spin component parallel to the field 
are pinned by H at low temperatures in a "one-in" and 
"two-out" spin configuration (see top right tetrahedron 
in Fig. nj. Because of the dipolar and exchange inter- 
actions, the spins on the a ^ a sublattices produce an 
internal local microscopic field, that enforces the 

two-in— two-out "spin ice rule" Q acting on the spins at 
sublattice a, hence forcing them to point "in" . A sim- 
ple estimate gives the (internal) Zeeman energy scale to 
be Hf£ = 2(5£> + Ji)/3 ~ 2.22 K at nearest neighbor 
distance, modified to 1.71 K using the Ewald method to 
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sum dipolar contributions [12|. In other words, the can- 
didate would-be field-decoupled z a = [111] spins, while 
experiencing zero applied external field H, are subject 
to an inward internal microscopic field of magnitude 
Wf{[ = (1.71K)/(7.09KT _1 ) = 0.24 T parallel to [111]. 

Such a large H™[ ~ 1.7 K precludes a phase transition 
at T <C iifii- However, a scenario for a low temperature 
phase transition can be reintroduced. Take the external 
applied field H canted away from perfect alignment per- 
pendicular to z a — [111] (see Fig. n H arrow), towards 
the inside of the tetrahedron (i.e. towards the [110] direc- 
tion), at an angle 9 such that 0(|H|) ss Hjj^/|H| in the 
limit of large applied field strength. For this orientation 
of H, the net local (internal + external) field acting on 
the a sublattice geometrically vanishes, so that the spins 
on a become (free) paramagnetic species. Spins on this 
a FCC sublattice are actually third nearest-neighbors 
on the pyrochlore structure (see Fig. 1), and therefore 
predominantly interact among themselves via long-range 
dipolar interactions (assuming momentarily zero third 
nearest- neighbor exchange J3). Interestingly, this field- 
decoupling scheme freezes out (decimates) spins on the 
three a ^ a sublattices at large field, and consequently 
generates a very rare physical realization of an FCC Ising 
(dipolar) magnet |2fij . 

We now determine the ground state for the spins 
on the a sublattice. First, recall that dipolar-coupled 
three-component (Heisenberg) spins have a ferromag- 
netic ground state on the FCC lattice, with an energy 
that is independent of the magnetization direction (as 
for all cubic lattices) |2l|. Consequently, with the field- 
decoupled a spins above residing on an FCC sublattice 
with Ising direction, z a = [111], we expect a phase tran- 
sition to a ferromagnetic ground state, with the bulk 
magnetization along ±[111] . Results from Monte Carlo 
simulations considering only spins on the a FCC sublat- 
tice, shown in Fig. |2 confirm that the ordered state is 
indeed ferromagnetic below 0.65 K. The inset of Fig. 
shows the results of simulations on the full pyrochlore 
lattice, with external fields ranging from 0.5 T to 1.5 T 
and missaligned from [112] with an angle 0(|H|), so as to 
give the geometrical cancellation above. There is no dis- 
cernable difference between the results above 1.5 T and 
the results obtained when considering only the a sub- 
lattice (labelled "Inf."), confirming the simple geometric 
arguments above. 

The transition we observe is at ~ 0.65 K, not 0.35 
K as found in powder specific heat measurements ||. 
However, adding a very small antiferromagnetic third 
nearest-neighbor exchange J3 = —0.022 K gives a tran- 
sition to a ferromagnetic ground state at a temperature 
close to the experimental value of 0.35 K. Interestingly, 
an antiferromagnetic J3 of this magnitude (— 0.03K < 
J3 < — 0.01K) appears able to explain the paramagnetic 
zone boundar y sc attering observed in neutron scattering 
of Dy 2 Ti 2 7 lUj. A more antiferromagnetic (negative) 




FIG. 3: Specific heat from MC simulations of (3 chains, equiv- 
alent to the limit large of [110] H applied to a single crystal 
of Dy2Ti207. The inset shows the q = X order parameter 
obtained via loop Monte Carlo (loop, open squares) [T^ . |23|| 
for a [110] field of 0.05 T and compared with conventional 
Monte Carlo results using single spin flip (SSF, filled circle). 

J3 changes the ordered state to a Neel state with ordering 
wavevector q = (1/2, 1/2, 0). The present discussion sug- 
gests that measurements on a well oriented single crystal 
could allow for a precise determination of a very small 
J3 value. A subsequent analysis of either DC suscepti- 
bility 0| , specific heat or elastic neutron scattering 
[l7)| could then help determine the full set of exchange 
parameters J± — J 2 — J3 in spin ice materials. 

The above calculation of assumed static spins on 
the three a ^ a sublattices, and neglected thermal fluc- 
tuations. Conversely, the applied "compensating" field 
parallel to [111] from the canted H(6>) is static. Hence, 
for any finite large [112] external field plus weak [111] 
compensating field of approximately 0.24 T, there is at 
nonzero temperature a net "overcompensating" outward 
field along [111] due to the thermal fluctuation (spin flips) 
of spins on the three a ^ a sublattices. The global (111) 
Ising symmetry is therefore always explicitely broken. 
The 9 — H relationship above is akin to the first order 
(H = 0, T < T c ) transition line in a conventional Ising 
ferromagnet. By correcting the angle 0(|H|), evaluated 
above for T = 0, to the one that would cancel the average 
Hjj£ at T = T c , a true spontaneous symmetry breaking 
can then occur at T c . For J3 < —0.03, the overcompen- 
sating field is not conjugate to the Neel order parameter 
and the transition is thermodynamically sharp for large 
applied H perpendicular to [111]. 

We now turn to the special case in which a second sub- 
lattice becomes accidentally field-decoupled as the field 
is rotated in the plane perpendicular to one of the four 
z a directions. This occurs for H along the [110 ] and re- 
lated crystallographic directions Q, @, E, E3i E-l E3 • 
As illustrated in Fig. ^ the [HO] field direction couples 
to the Ising spins on two FCC sublattices (forming the 
a chains), but leaves spins on the two FCC sublattices 
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forming the (3 chains fully decoupled. Unlike the generic 
case above with only one (a) decoupled sublattice, the net 
local internal field from static spins on the a chains van- 
ishes on lattice sites that belong to (3 chains. Although 
the spins on the /3 chains are then in principle free to 
undergo an ordering transition via some long 
change or dipolar interaction, experimental 
and numerical 0, evidence for this transition has 
so far been elusive. Our Monte Carlo simulation results 
on the dipolar spin ice model show strongevidence for 
a transition to the so-called q = X state 1] for H par- 
allel to [110], where ferromagnetic one-dimensional (ID) 
(3 chains along [110] orient antiferromagntically with re- 
spect to their four nearest neighbour chains (see Fig- EJ - 
This is found by studying both the q = X order param- 
eter 0] an( i relative zero-temperature energies of the 
various competing spin ice ground states 0, Q. A 
transition to q = X occurs for all fields greater than 
- 0.02 T, although for weak fields (H < 0.10 T) non- 
local loop spin flips are needed in the simulation to find 
the ground state jl^]. For fields weaker than 0.02 T, 
the transition is to the q = (0, 0, 1) long-range ordered 
state found by Melko et al. |l2|, |23j . The critical tem- 
perature T c 110 (|Hno|) for the q = X ordering is strongly 
field dependent for field strengths between 0.02 T and 
0.2 T [HEl, 

but saturates at a critical temperature of 
T c 110 (oo) ~ 1.2 K for H > 0.2 T, as also recently found in 
Ref. Jl4| , but contrary to what was obtained in Ref . . 
For large H along [110] the transition to q =X proceeds 
via a ID to 3D crossover with a correlation length £i£> 
of approximately 20 nearest-neighbor lattice spacing at 
the crossover temperature. This crossover at such a large 
£id makes the study of finite-size scaling and 3D critical 
behavior very demanding at this transition. The pres- 
ence of dilute defects or impurities or the lack of perfect 
alignment 0] of H along the [110] direction might be 
at play in preventing the development of perfect q = X 
long range order in experiments [l|, [l3|, |22| • Finally, we 
note that introducing J3 = —0.022 K in Eq. (1) results 
in a shift of T 110 (oc) to T c 110 (oo) - 1.6 K, in disagree- 
ment with the position of the specific heat feature at 1.2 
K Such a large shift of T c 110 is consistent with 

the T c 110 calculated using a ID chain mean-field theory. 
Overall, the results presented in this paper suggest that 
there might be important competition between J\ , J2 and 
J3 exchange in Dy2Ti20y, and that the original descrip- 
tion of D y 2 Ti2 O 7 using only a nearest-neighbor exchange 
Ji = -3.72 K [lllnl should be revisited. 

In conclusion, we have identified a mechanism that 
causes field-independent phase transitions in the dipolar 
spin ice model when subject to a large external applied 
field. For polycrystalline Dy2Ti207, we suggest that the 



transition at J 1 ] 112 ' ss 0.35 K can arise from crystallites 
for which the field is accidentally quasi-perpendicular to 
one of the (111) directions, but slightly canted towards 
the corresponding (111) direction to cancel internal ex- 
change plus dipolar local field. This proposed mechanism 
could be verified in a single crystal experiment. We also 
postulate a small third neighbor antiferromagnetic ex- 
change to account for the precise tJ 112 '. In addition, we 
found unequivocal support for the existence of a tran- 
sition to a long-range ordered q = X state for dipolar 
spin ice in a magnetic field along [110] in excess of 0.02 
T, which is postulated to account for the tJ 110 ' 1.2K 
transition in polycrystalline Dy2Ti207. As an interest- 
ing by-product of this work, we have identified a novel 
realization of an Ising FCC dipolar ferromagnet. 
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